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ABSTRACT 

The European Space Agency's Gaia mission is scheduled for launch in 2013. It will 
operate at L2 for 5 years, rotating slowly to scan the sky so that its two optical 
telescopes will repeatedly observe more than one billion stars. The resulting data 
set will be iteratively reduced to solve for the position, parallax and proper motion of 
every observed star. The focal plane contains 106 large area silicon CCDs continuously 
operating in a mode where the line transfer rate and the satellite rotation are in 
synchronisation. 

One of the greatest challenges facing the mission is radiation damage to the CCDs 
which will cause charge deferral and image shape distortion. This is particularly im- 
portant because of the extreme accuracy requirements of the mission. Despite steps 
taken at hardware level to minimise the effects of radiation, the residual distortion 
will need to be calibrated during the pipeline data processing. Due to the volume and 
inhomogeneity of data involved, this requires a model which describes the effects of 
the radiation damage which is physically realistic, yet fast enough to implement in 
the pipeline. The resulting charge distortion model was developed specifically for the 
Gaia CCD operating mode. However, a generalised version is presented in this paper 
and this has already been applied in a broader context, for example to investigate the 
impact of radiation damage on the Euclid dark-energy mission data. 
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; 1 INTRODUCTION 

Gaia is a European Space Agency mission, planned for 
launch in 2013, which aims to perform global astrometry on 
approximately one percent of the estimate d galactic stellar 
population to unprecedente d accurac^F]f see [Perrvman et al.l 
I2OOII : iLindegren et~al]|2008l ) . The scanning satellite will op- 
erate at the earth/moon-sun Lagrangian point for 5 years, 
rotating slowly so that its two optical telescopes will re- 
peatedly observe more than a thousand million stars. The 
resulting dataset will be iteratively reduced to solve for the 
position, parallax and proper motion of every observed star 
(additionally, dispersed spectra will be obtained for a subset 
of these objects). The mission requirements for the astromet- 
ric precision are in the micro- arcsecond ([J.as) regime. These 
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^ For science performance predictions see Ide Bruiinj 

1I2OI2I ). for current updates on the science performances see 
|http:/ /www.rssd.esa.int/gaia| 
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stringent astrometric requirements translate into milli-pixel- 
level accuracies on the stellar image location per individual 
CCD observation, e.g., [B-od'hommc_tit_aL., (201 1) . 

The Gaia focal plane (e.g.. Ide Bruiine et al.|[2010l h the 
largest to date to be sent to space, contains 106 large- 
area silicon CCDs continuously operating in Time Delay 
Integration (TDI) mode at a line-clocking rate which is 
in synchronisation with the satellite rotation rate. Each 
CCD was custom-made by e2v technologies for the Gaia 
project and consists of 4500 x 1966 (parallel x serial) pix- 
els of physical dimension 10 x 30 pm (for further detail see 
iProd'homme et al.|[20Tll ). Due to a combination of the large 
number of transfers of the charge packet (parallel) and fast 
clocking speed (serial), the limits placed on the amount of 
radiation shielding due to mass-budget constraints, the pre- 
dicted radiation environment at L2, and the extreme accu- 
racy requirements on the image location estimates. Charge 
Transfer Inefficiency (CTI) of the CCDs was identified as a 
challenge for the Gaia mission at an early stage. 

CTI during the transfer of electron packets from pixel 
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Figure 1. Upper Panel: Simulation of a 5 X oversampled image 
in the abscence of CTI. Lower Panel: The effect of parallel CTI 
on the image in the top panel. The CTI effect was simulated 
using the model described in this paper. Note that unrealistic 
model parameters were chosen in order to exaggerate the effect 
for clarity. The readout direction is to the right. 



to pixel is due to the stochastic capture and release of signal 
electrons into and out of trapping sites in the silicon lattice 
structure. This will have two major effects upon Gaia data: 

(i) the removal of charge from the signal within the win- 
dowed imagf0 resulting in a loss of signal and a correspond- 
ing irretrievable degradation in the signal-to-noise ratio. 

(ii) A distortion and centroid shift of the image due to the 
emission of electrons into charge packets different to those 
which they were captured from. 

These effects are schematically illustrated in Figure [1] 
and Figure [2l 

The distortion of the image shape, if uncalibrated, will 
introduce a systematic bias on the measured image loca- 
tion and result in a degradation of the end-of-mission as- 
trometric accur acies from [xas to mas leve ls. Indeed, based 
on simulations, iProd 'homme et al.l (12011') estimate a bias 
of up to ~ 0.2 pixels in the CCD location measurement 
(~ 10 mill-arcseconds) as a result of CTI in the uncali- 
brated case, which is in agreement with on-ground test data. 



^ In order to meet telemetry constraints and to minimise readout 
noise, windows are placed around each confirmed source and only 
data from within these windows arc read out and sent to ground. 
The window dimensions depend on the magnitude of the source. 
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Figure 2. Simulation of a 5 X oversampled image that is binned 
in the serial direction (the vast majority of Gaia observations will 
have on-chip binning applied in the serial direction). The effect 
of parallel CTI on the image is plotted in dashed. The CTI effect 
was simulated using the model described in this paper. Note that 
unrealistic model parameters were chosen in order to exaggerate 
the effect for clarity. The readout direction is to the right. 



The main reasons that the effect on the performance is so 
large is due to the harsh radiation environment at L2 (dom- 
inated by energetic solar protons) and the stringent mis- 
sion requirements on the image location measurements. It is 
well established that the interaction of high-energy protons 
with CCD devices can result in displacement damage in the 
silicon structure and the generation of 'new' energy levels 
in the semiconductor band-gap, otherwise known as 'traps', 
that can capture and, after some time, release the electrons. 
These radiation-induced traps dramatically increase the CTI 
of the detector and raise the radiation issue to being mission- 
critical in the absence of adequate mi tigation. For furthe r 
discussion of radiation-induced CTI see iHardv et al] (|l998D : 
IJanesickl (|200ll V 

In order to minimise the CTI effects, a number of hard- 
ware solutions has been implemented for Gaia. For example, 
a mechanism for the periodic injection of charge into lines of 
the CCD in order to keep traps filled has been implementecffl 
(charge injection). A supplementary buried channel has also 
been incorporated into the CCD design to funnel charge 
from small electron packets into a smaller pixel volume dur- 
ing transfer in order to minimise the number of traps that 
will be encountered at low signals (for further reading see 
ISeabroke eralll2013l ). However, in order for Gaia to meet its 
scientific requirements, a calibration of the residual effects 
of radiation-induced CTI is mandatory. 

The entire set of Gaia observations will be used to solve 
for the astrometric parameters of each star, as well as the 
satellite attitude parameters and calibration parameters us- 
ing the Astrometric Global Iterative Solution (AGIS) de- 



Note that these periodic charge injections will also act to de- 
crease the scientific quality of the data in two main ways. Some 
observations will be corrupted by the injection of charge across 
the acquired windowed data. Also, the trail from the injection 
will provide an additional background and thus inject additional 
noise into the data. However, the benefits of the implementation 
of periodic c harge injection outweigh these downsides (see t he 
discussions in lProd'homme et al.l l|20l3) and lHoU et"aLl 1I2OI2I) '). 
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scribed in iLindegren et alj (|2012l '). Radiation-induced distor- 
tions of the images are a source of systematic error which, 
given a suitable model, could be removed as part of the 
global solution. Therefore, the current baseline approach to 
tackle the calibration of radiation-induced trapping effects in 
Gaia data is through the application of a forward modelling 
technique in the pipeline processing. Under this scheme, the 
parameters of the CTI-free image (location and integrated 
flux) and the parameters of a model which can replicate 
the CTI effects on the CTI-free image are iteratively ad- 
justed until the predicted distorted image best matches the 
observed image. This approach requires a model that can 
replicate the distortion of a sampled image due to radiation- 
induced trapping effects, we call this a Charge Distortion 
Model (CDM). The forward-modelling approach is prefer- 
able to a corrective approach in order to preserve the noise 
properties of the data, as well as having the flexibility to 
treat dispersed spectra and non-point source objects in the 
same manner as isolated single stars in the data processing. 
However, the CDM must be realistic and flexible enough to 
reproduce the distortion effectively, whilst simultaneously 
being fast enough for implementation in the data processing 
pipeline. The CDM can be empirical and need not be based 
upon detailed physical modelling of electron trapping and 
release. However, any model based upon physical consider- 
ations is certain to be superior in terms of application over 
a broad parameter space using the fewest possible variables. 
Since computational speed is of high concern, the model 
need not treat individual pixel transfers, but instead, ap- 
proximate the trapping and release during the transit of an 
entire CCD column (or, indeed, the readout register) in a 
single calculation. In order to treat the non-integrating (non- 
TDI) charge injection lines, transfer through the serial regis- 
ter and also for general applicability, the model should also 
have the capability to treat the charge transfer when the 
signal is non-integrating (imaging mode). 

During the course of Gaia development and preparation 
studies, a great deal has been learned regarding the phe- 
nomenological interpretation of trapping, de-trapping and 
electron confinement from the very large to the very small 
signal regimes. Over the course of five radiation test cam- 
paigns, many datasets have been acquired in TDI mode us- 
ing partially-irradiated Gaia CCD s in order to characterise 
the CTI effects. Monte-Carlo fe.g.. |Prod'homme et alll201ll ') 
and analytical models have also been developed which go 
some way to reproducing Gaia radiation test data. However, 
the sub-pixel level Monte Carlo CCD models developed for 
Gaia are far too computationally expensive to be applied 
as a CDM in the forward modelling scheme which needs 
to run on ~ 10^^ separate CCD observations collected over 
the course of the missiorQ In this paper, we present a model, 
based upon a physical understanding of trapping effects, but 
with some very significant simplifications such that it can be 



* In order to demonstrate the speed advantages of this model, 
a typical Gaia CCD observation (6 X 12 pixel image) was run 
through a Java implementation of this CDM (not optimised for 
speed). With 7 trap species used by the model the it takes ~ 0.25 
milli-seconds to process using a standard laptop with a Intel Core 
i5 processing chip. Running the sam e image through a Monte 
Carlo model jProd'homme et al.ll201lh takes 60 seconds. 



incorporated into the data processing pipeline. Further dis- 
cussions and quantitative predictions on the performance of 
Gaia in the presence of CTI (including t he use of this CDM 
for CT I-mi tigation) can be f ound within [Prod'homme et ah! 
(|201ll ) and lHoll et al] l|2012l ). 

The paper is structured as follows. In Section [2] we ad- 
dress the constraints imposed on the model from the analysis 
of Gaia-specific test data acquired in TDI mode. The details 
of the model are then presented in Section [S] In Section |4] 
we show an overview of the application of the model to test 
data and, finally, in Section [5] we present our conclusions. 



2 BACKGROUND TO THE MODEL 

The objective is to derive an analytical CDM that takes into 
account recent illumination and charge injection historjjf] 
and is based upon acce pted physical theo ry, i.e. Shockley- 
Read-Hall theory (e.g.. [Hardv et al]|l998l ) and knowledge 
gained from Gaia testing. The first iteration of such a CDM 
was proposed by iLindegrenI (|2008l ). however, his model is 
not based on physical theory and was always intended as a 
place-holder. 

In order to apply the Shockley-Read-Hall equations it 
is necessary to model the volume and density of the elec- 
tron packets as they move through the CCD. This will give 
the numbers of traps encountered by an electron packet ac- 
cording to the trap density, as well as the probability that 
a trap captures or releases an electron. In this respect, the 
models developed for Gaia fall into two broad categories: 
Gonfinement Volume (or volume driven) models and Den- 
sity Distribution (or density driven) models. 

CCD models usually consider electrons to be contained 
in charge packets resting in potential minima defined by the 
applied electrode potentials. The electron cloud is typically 
assumed to have a finite confinement volume which increases 
as more electrons are added. Hence a model such as: 



Vg \FWCJ ' ' 

might be adopted for the volume of the charge cloud 
(Vc), where Nc, is the number of electrons in a pixel, FWC 
is the pixel Full Well Capacity in electrons and Vg is the 
assumed maximum geometrical volume that electrons can 
occupy within a pixel (i.e. the volume of the electron cloud 
when = FWCfl In most models, the value of /? would 



^ The number of electron-trapping sites that can capture an elec- 
tron will of course depend upon the fraction of these sites that 
are currently vacant, the trap occupancy. The trap occupancy will 
thus depend on factors such as the time since last charge injec- 
tion, the constant diffuse optical background incident upon the 
detector, the recent transit history of other light sources across 
the pixels in question etc. 

^ A different function describing the behaviour of the confinement 
volume with electron packet size could also be used which could be 
optimised for a specific CCD pixel architecture, however here we 
present a generalised model. In fact, this function could possibly 
be optimised for use in the Gaia data processing after analysis of 
real CTI-affected Gaia data at some stage during the mission. 
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be assumed to be close to 1, the density of confined elec- 
trons would be almost constant and high enough for trap- 
ping within the electron packet volume to be considered in- 
stantaneous. We may refer to models with these properties 
as volume driven models where all vacant traps within the 
electron packet will capture an electron whilst those outside 
cannot, and the amount of trapping observed will be driven 
by the electron packet volume. Early models developed for 
Gaia were all volume driven. However, there was a major 
shortcoming of this assumption which causes them to fail 
for Gaia. 

Volume driven models are unable to reproduce or to ex- 
plain the observed effect that a small number of background 
electrons can have upon overall trap occupancy. It is ob- 
served that very low level optical background (fat zero) is 
able to suppress a disproportiona t e fraction o f slow traps 
(e.g.. iBrown fc van LeeuwenI l2008l : IShortll2006h . According 
to volume driven models, this would require that background 
electrons must occupy a volume many thousands of times 
larger than an equivalent number of signal electrons. In ad- 
dition, Gaia testing clearly indicates that a diffuse optical 
background preferentially eliminates traps with long release 
time constants (slow traps) but does little or nothing to 
suppress faster traps. The effect of background upon the 
traps suggests that a different model must be applied allow- 
ing electrons to reach the traps in the first place. In other 
words, we must assume that jS in equation [T] is closer to 
rather than 1 and that it is primarily the density of electrons 
rather than the volume which increases as more electrons are 
added. In this way, all electron packets will encounter a sim- 
ilar population of traps but the probability that a given trap 
captures an electron in a given time will depend upon the 
local electron density which increases with signal size. 



3 MODEL DESCRIPTION 

In this section we break down the model description into four 
differerent sections. We describe the modelling of the elec- 
tron trapping in TDI mode, the trapping in imaging mode, 
electon release and, finally, the tracking of the trap occu- 
pancies as a function of tim^lJ. 



3.1 Electron trapping in TDI mode 

3.1.1 Interaction volume 

We take equation [1] as a starting point, noting that /3 is not 
constrained to values close to 1 but may take any value be- 
tween and 1 allowing complete flexibility with regard to 
the question of volume or density driven models. Clearly, if 




Figure 3. Variation of electron confinement volume, Vc with 
number of electrons, A^e according to equation [T] 



is entirely driven by the interaction volume of the electrons 
with the traps. Figure |3] illustrates the variation of electron 
confinement volume, Vc as a function of A''e and /3 (according 
to equation [T|. 

Since we intend to treat the trapping along the column 
in a single step, we may derive the fraction of the total vol- 
ume of silicon in the column that was "seen" by a given 
sample (5*), by integrating equation [T] as the signal electrons 
accumulate along the CCD columijf] from TVe = to A'^e = 5": 



volume that sees electrons 
volume of column 



\FWC 



1 

1 

TT/3 \FWC 



ANc 



S 



(2) 



where it should be understood that S, the input to the 
CDM is the size of the sample (in electrons) that would be 
read out of the CCD in the absence of trapping, i.e. the un- 
distorted sample. Assuming that Nt traps are distributed 
uniformly throughout the volume of the column, then the 
fraction of these traps with which the 5* electrons can inter- 
act is equal to the interacting fraction of the column volume. 
However, if we assume that No of the traps within the col- 
umn are already occupied, then we can write the fraction of 
all traps in the column that will interact with the S electrons 
and that are vacanl|f|: 



Ft 



vacant traps that see electrons 
total number of traps 

1 ( s y No 

I + 13 \FWC J Nt 



(3) 



tron density increases as electrons are added and trapping 
is entirely density driven. If P is equal to 1 then Vc grows 
linearly with A^e, electron density is constant and trapping 



^ It should be pointed out that since the model is analytical, any 
trap parameters derived with this model will not necessarily cor- 
respond to exact physical quantities, even when the model can 
acurately reproduce what is observed in the data. They should 
therefore be thought of as effective, rather than physical, param- 
eters. 



* By 'column' we mean the column of pixels in the parallel di- 
rection that an electron packet must be transferred along before 
reaching the readout register (also called the serial register). 
^ It should be noted that No describes the state of the traps in a 
column as seen by a given TDI sample and so refers to the state 
of each trap along the column at the moment in time when it is 
encountered by the current sample, i.e., immediately after being 
left by the previous sample. Therefore, it is not referring to the 
state of all of the traps in the column at the moment when a new 
sample begins to transit the CCD. 
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Similarly, we may write the fraction of signal electrons, 
S, that see vacant (rather than occupied) traps: 



electrons that see vacant traps 
total number of electrons 



1 



' _S_\l^ _ No_ 

> FWC J Nt 



1 

l+,8 



( s \t 
\FWC J 



(4) 



If trapping were allowed to proceed until all of the va- 
cant traps were occupied or all of the available electrons were 
captured, then the number of captures, Nc would be equal 
to the number of available traps or the number of available 
electrons (whichever is fewer). This may be expressed: 



N, 



c £— J-oo — 



NtFt + SFt 



1 

1+/9 



FWC J 



Xrhi 



FWC ) 



1 ( S \P 
1 + 3 V FWC I 



Nt 



1 

1+13 



No 



+ S 



( 1 ( s \P 

\l + fi\FWCl 

1+3 \ FWC j 



where Tc is the capture time constant: 

1 



(7) 



(8) 



The quantity a is the trap capture cross-section, is 
the electron density in the vicinity of the trap and vt is the 
electron thermal velocity: 



1 » 2 3 



Vt = 



mi 



(9) 



where is the effective electron mass in silicon which 
approximately equals half the free electron rest masj"l. 

If we assume that a number of electrons (A^e) are con- 
tained within an effective confinement volume (Vc), then: 



(JVtNe 

Substituting into equation [T] 



(10) 



-yS^ - No 



where 



Nt 



(1 + /3) FWCf 



(5) 



When implementing the model, it is more general to 
use a trap density, nt, rather than a number of traps in the 
column, Nt- Thus, substituting Nt ~ 2ntVgX where x is 
the number of TDI transfers or the column length in pixels, 
we get: 



2ntVgX 



(1 + P) FWCP 



(6) 



3.1.2 Finite trapping time 

Equation [5] may be used to model electron trapping in TDI 
mode in the event that capture times are always short in 
relation to electron-trap interaction times (i.e. semi instan- 
taneous trapping). However, it is well established through 
Gaia testing, that the very low electron densities being 
transferred in TDI mode give rise to electron capture times 
that are on the order of the pixel dwell time and it is there- 
fore not appropriate to consider instantaneous trapping. In- 
stead, we must evaluate the probability of electron capture 
according to electron density. In TDI mode the electron den- 
sity is increasing as the signal integrates along the CCD 
column. Therefore, most Gaia models evaluate the proba- 
bility of capture (and release) per line transfer or intra-pixel 
CCD phase (4 500 transfers or 18 000 phases, since each pixel 
contains four individual gate electrodes). However, for the 
purposes of this high speed model, we require an expres- 
sion giving the column-averaged probability of capture for 
a transit as a function of (integrating) signal size and trap 
occupancy. 

If we assume that traps in the CCD column interact 
with the electrons of a given sample for a period equal to 
half the TDI period {t) then the probability (Pc) that a 
vacant trap will capture an electron is: 



Pc = 1 — exp 



tavtNe 
2Vc 



and substituting for Vc from equation [T] 



Pc = 1 — exp 



tavtNc ( FWC V 
2Vg \ Nc ) 



(11) 



(12) 



Assuming that the TDI period, t, is constant, we then 



find: 



p,., , ^ -cNl-P , tavtFWC^ 
FcyJyie) = 1 — e ^ where a = — — 

(13) 

Equation [13] gives the capture probability (per vacant 
trap) as a function of the number of sample electrons, Nc. 
However, we require an effective capture probability for the 
entire column in which the number of electrons is increasing 
from Nc = to Nc — S where S is the size of the sample in 
the absence of trapping. To first order, we may evaluate an 
effective or column-average capture probability as a function 
of the input signal, S, using: 



column average PdNc) — Pc{S) 



- / Pc{Nc)dNc 
■5 Jo 



Sjo 



l-e-"''^' ^ dNc (14) 



which yields: 



Particles in a specific solid material, and under the influence of 
an external electromagnetic field, can generally be approximated 
to behave as if they were free particles under semi-classical solid- 
state physics, however, a modified mass needs to be used, which 
is the effective mass for that specific material. 



© RAS, MNRAS 000,[T]-?? 



6 A. Short et al. 





P=0.9 
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Figure 4. Column-average electron capture probability, Pc{S), 
according to equation 1151 (black) and the simple approximation 
given in equation 1181 (red) . Note that a rather small value of cap- 
ture cross-section (3 X 10~^* cm^) has been used in this example 
to illustrate the dependence over a large signal range. Note that 
the difference is smaller for larger values of the capture cross sec- 
tion where the curve would be almost flat and close to unity. 



-yS^ ^ No 



1 — CXp < — Q 



1-/3 



(19) 



Note that equation [19] can give negative trapping when 
the number of occupied traps exceeds the number of inter- 
acting electrons and vacant traps. If the de-trapping process 
were the reverse of the trapping process, then this could be 
used directly to model trapping and de-trapping. However, 
trapping and de-trapping time constants are very different 
so this is not the case. When implementing equation 1191 it 
is therefore necessary to avoid negative trapping. Thus, in 
summary, we have a model for electron capture by one trap 
species in TDI mode: 



Nc = MAX 



where 7 : 



1-/3 



2ntVgX 



(1 + 13) FWCf' 



and a = 



, 
(20) 

tavtFWC^ 



2Va 



(21) 



1 - 



-r(i + ^) (as-0^^")r(i. 



13-1 ' "'^ 



where 

r(^): 

is the gamma function and 



' du 



(15) 



(16) 



(17) 



is the incomplete gamma function. Equation 1151 may be 
used in a model but it is unwieldy and slow computation- 
ally to implement, especially when attempting to evaluate 
the initial (or equilibrium) trap occupancy prior to the tran- 
sit of the first signal paclcet. Instead, we make a very simple 
approximation that Pc{S) « Pc{S/2) and then test the va- 
lidity by comparison with equation [151 Thus 



assume Pc(S) ~ 1 — exp 



(18) 



Figure [4] illustrates equation [15] (black) and the simple 
approximation given in equation 1181 (red). The agreement 
is certainly good enough to justify the use of the simpler 
expression. Hence, we may combine equations [5] and [18] to 
give an expression for the number of electrons, Nc, that will 
be captured by a given trap species during the transit of an 
integrating signal packet, S, as a function of physical trap 
parameters and the initial trap occupancy. No'. 



3.2 Electron trapping in imaging mode (non-TDI) 

There are several circumstances where it is necessary to 
model trapping in non-TDI mode. These include the treat- 
ment of charge-injection lines which are injected at the top 
of the CCD and do not integrate along the CCD, modelling 
transfer along the serial register and of course applying the 
model to other missions which are not operating in TDI 
mode at all. 

There are only a couple of differences in the final equa- 
tions related to the fact that the signal is not integrating 
during the transfer. In imaging mode, the model may be 
summarised: 



Nc = MAX 



jS^ - N o 



L(l_exp{-a5^-^}) 



(22) 



where 



7 -- 



2ntVg 
FWCP 



[x + i) 



and 



tavtFWC^ 
2K 



(23) 

where x indicates the position on the CCD (in pixels) 
of the image to be processed, i.e., small x for an image close 
to the readout node and large x for objects far from the 
readout node and i is the pixel coordinate of each sample 
within the image to be processed. 



3.3 Electron release 

If we assume that every occupied trap in the column has 
a time equal to one TDI period it) in which to release the 
captured electron into a given sample, then the probability, 
Pr, that the trap will release the electron into the sample is 
simply: 



Pr = 1 



(24) 
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where Tr is the trap release time constant. The number 
of electrons released into the sample during a transit along 
the column is then simply: 



iV, = iVo|l-e-^| 



(25) 



1 — exp 



A = 



1-/3 



C = l-e~ 



75f + 1 
and 



and 



(28) 



3.4 Column state parameter and initial state 

In implementing any CCD trapping model, it is very im- 
portant to consider the state or occupancy of the traps (the 
value of No) at the beginning of the simulation. In the event 
that a continuous, unbroken stream of realistic data (in- 
cluding background etc.) were to be simulated (much longer 
than the capture or release time constants of the slowest 
traps), it would be sufficient to start with all traps empty 
or all traps full and to disregard the first few minutes of 
simulated data rather than to explicitly initialize No with a 
realistic value. However, this is rarely practical, the state of 
the slower traps prior to the transit of a source will be domi- 
nated by two main factors. These are the time since the last 
charge injection and the level of constant or slowly varying 
optical background (some fraction of slow traps kept full). 
We may therefore treat these two dependencies explicitly in 
order to quickly determine the state of traps prior to the 
transit of an astronomical source (or test data image). 



3.4-1 Diffuse optical background 

The effect of a constant (or slowly varying) low-level opti- 
cal background upon the occupancy of slow traps has been 
demonstrated during the Gaia radiation test campaigns (see 
Section [2]). This reveals information about the temporal de- 
pendency or dynamics of trapping since the exposure time 
of traps to a constant background is radically different from 
the exposure time to transient sources. It also reveals some- 
thing about the confinement volume of electrons within pix- 
els since low levels of background appear capable of filling 
entire trap species. 

Considering TDI mode for Gaia, the initial trap oc- 
cupancy, No\init due to a constant optical background is 
determined as follows. It is assumed that in the continuous 
presence of background electrons, traps will reach an equi- 
librium state of occupancy when Nc — Nr. Thus: 




where Sdob is the diffuse optical background level in 
electrons per pixel measured at the readout register. 
Solving equation 1261 for No\ init gives: 



3.4.2 Time since charge injection 

In the case of a high level of charge injection it may be suf- 
ficient to assume that the injection fills all traps which then 
subsequently empty at a rate determined by the exponential 
trap release time constant until they reach the occupancy 
level determined by the diffuse optical background. For ex- 
ample, adding an exponential term to eQuation l27l the initial 
state of the column would be given by: 

NoUmt^2ntVgXe-— + ^^{l~e- — j (29) 

where td is the time elapsed since charge injection and 
Tr is the trap release time constant. 

However, in the case of lower levels of charge injection 
where it may not be adequate to assume that all traps are 
saturated, it is better to simulate the capture of electrons 
from the charge injection lines explicitly to determine the 
fraction of traps that are filled. These then empty according 
to the exponential trap release time constant as before. In 
order to simulate charge injection correctly, it is necessary 
to consider that charge injection is effectively in image mode 
rather than TDI mode. 



4 EXAMPLES OF MODELLING RESULTS 

It is beyond the scope of this work to provide a detailed 
analysis of the application of the model to Gaia test data. 
Instead, we provide a limited graphical indication that the 
model is qualitatively capable of reproducing trapping ef- 
fects observed in test data. 

Figure [5] (left panel) indicates the distortion observed 
in the parallel direction in Gaia test data as a function of 
time since charge injection. The datapoints are measured 
data and the solid lines are CDM model results. Note that a 
single set of trap parameter values was used to generate all 
of these model curves. Figure [S] (right panel) shows the fit 
of the model to radiation damaged dispersed spectral data 
(similar to what will be observed with the Radial Veloc- 
ity Spectrometer (RVS) instrument onboard Gaia). Curves 
are shown for data taken on non-irradiated and irradiated 
sections of a CCD, with the output of the CDM also over- 
plotted. Again, a single set of model parameters has been 
selected giving a reasonable fit to all data. The CDM has 
also been used to exami ne the effects of radia t ion damage on 
CTI in Euclid (e.g., see lLaureiis et al]|2012l : ICropper et al.l 
l2013h CCDs. 



No\ 



AB 

A + C 



(27) 



where: 



5 CONCLUSIONS 

The ESA astrometric mission Gaia, scheduled for launch in 
2013, will be subject to a systematic bias in its measure- 



© RAS, MNRAS 000,[T]-?? 



8 A. Short et al. 



ISOOO 







Time since charge 
injection 






-non-irradiated 






1000 TD! periods 
-model 


'It / 

■ 1 / 


w • 

\\ 


5000 TDI periods 
-model 


Hi 

■■11 / 




30000 TDI periods 


In 




- model 


ll 1/ 


\ Vi • 


120000 TDI periods 
-model 



300,0 



- non-irradiated data 



-irradiated data 




-5 -4 -3 -2 -1 1 2 3 4 5 
TDI pixels 

Figure 5. Left: Rcprcscntivc fits to Gaia test data for Gaia magnitude ~15 point sources with varying time since cliarge injection. The 
non-irradiated curve is derived from data acquired on a section of a CCD which has not undergone radiation damage. The other curves 
are derived from data acquired from a radiation-damaged section of the same detector, with a uniform damage level throughout this 
section. Right: Similarly, a section of Gaia dispersed test data for Gaia magnitude ~10 and with a diffuse optical background of 0.77 
e~pix~^s~^ . These data are representative of what will be observed with the Radial Velocity Spectrometer (RVS) instrument onboard 
Gaia. For the fits shown in both plots model parameters are derived using a global fit to all test data from magnitude 13.25 to magnitude 
20 and with delays since prior charge-injections from 30 ms to 120 s. Note that readout direction is to the left for both plots. 



ments due to radiation damage-induced CTI in its CCD 
detectors. In order to sufficiently calibrate out this bias it 
is planned to forward model the charge distortion to the 
CCD images at the CCD sample level. The use of an an- 
alytical charge distortion model, which must be computa- 
tionally inexpensive, is fundamental to Gaia reaching its 
scientific goals. It treats electron trapping and de-trapping 
within each CCD column (and serial register) in a single 
step ensuring that it is computationally fast, although still 
based on physical theory. A description of the model is pre- 
sented, and its ability to qualitatively reproduce the CTI 
effects in Gaia test data is demonstrated. The model can be 
applied to CCDs operating in both TDI and imaging mode 
and hence is fully general. 
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